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ABSTRACT. It is shown that the vacuum Einstein equations for an arbitrary stationary axisym- 
metric space-time can be completely separated by re-formulating the Ernst equation and its associated 
linear system in terms of a non-autonomous Schlesinger-type dynamical system. The conformal fac- 
tor of the metric coincides (up to some explicitly computable factor) with the r-function of the Ernst 
equation in the presence of finitely many regular singularities. We also present a canonical formulation 
of these results, which is based on a "two-time" Hamiltonian approach, and which opens new avenues 
for the quantization of such systems. 



Introduction. In this letter we demonstrate that the vacuum Einstein equations for space-times 
with two commuting Killing vectors can be re-formulated in terms of a pair of compatible ordinary 
matrix differential equations. Similar results can be shown to hold for the more general equations 
obtained by dimensional reduction from higher-dimensional theories of gravity and supergravity with 
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matter couplings to two dimensions. As a by-product, we establish a previously unknown relation 
between the conformal factor of the associated metric and the so-called r-function, which plays a 
pivotal role in the modern formulation of integrable systems ^. Thirdly, we present a canonical 
formulation of these results, which avoids certain technical difficulties encountered in previous treat- 
ments. Our results suggest that an exact quantization of axisymmetric stationary (matter-coupled) 
gravity by exploiting techniques developed for flat space (quantum) integrable systems j^, Q is now 
within reach. 

The Ernst equation and related linear system. We start from the following metric on 
stationary axisymmetric space-time 

ds"^ = f-^e^^idx"^ + dp^) + p^dcj)'^] - f{dt + Fd<i)f (1) 



where (x, p) are Weyl canonical coordinates; t and (j) are the time and angular coordinates, respectively. 
The functions f{x,p), F{x,p) and k{x,p) are related to the complex Ernst potential £{x,p) by 



f = Re£ , F^ = 2p ^l ^ f l , k^ = 2ip 



{£ + £f 



{£ + £?' 



(2) 



where ^ = x + ip , = x — ip; hereafter subscripts ^, ^ denote partial derivatives with respect to these 
variables. In terms of the potential £{C,C)j Einstein's equations for the metric (||) in particular imply 
the Ernst equation |^] 

((^ - Og^g-')^ + ((^ - Og^g-')^ = o (3) 



with the symmetric matrix 



£ + £ 



^ 2 i{£-£)^ 



^ i{£ - £) 2££ ^ 



(4) 
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The equation (^) for k{^,^) may be equivalently written in the form 



d log h = q (5) 
where h := e^^ is the conformal factor. Using (^) one can show that the one- form q defined by 

q = ^tr(<755-')'dC + ^t^[g-^g-^fdi (6) 

is closed, i.e. dq = 0. Equation @ is the compatibility condition of the following linear system |^]: 



f = ^* «=y* (7) 



where 

U='-^ y = ^; (8) 
1-7 1+7 

and ^'(7, ^, ^) is a 2 X 2 matrix, from which the Ernst potential and thus the metric (|l|) can be recon- 
structed. The function 7(^, ^) is a "variable spectral parameter" subject to the following (compatible) 
first order equations 

7 _1 +7 __ _7 1-7 .gX 

They are solved by 



7(C,C» = 1^ - ^ - (10) 

with If G C a constant of integration, which can be regarded as the "hidden" constant spectral 
parameter. For the linear system (|^, we can use either 7 or when 7 is expressed as a function of 
w according to (|^), the linear system (^) lives on the two-sheeted Riemann surface of the function 



{w — Cj{w — S^). Both the constant and the variable spectral parameters w and 7 are needed for a 
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proper understanding of the (infinite-dimensional) liidden symmetries of (^) |^ |9|. Furthermore, 



d _ d 7^1 + 7 ^ _ g_ 7 1 - 7 ^ 



The poles of (|^ in the complex 7-plane are thus produced by differentiation of 7 according to @. 
Choosing {'j,^,^) as independent variables, we get the following relations from (0) and |10|: 



9^9 ^ = -^^1^ ^\i=i 9^9 ^ = -^^1^ ^l7=-i (12) 



where the subscript 7 denotes differentiation with respect to 7. 

Deformation equations. Although we shall keep in mind the Ernst equation and its associated 
linear system @ as our principal example, the results to described below hold for arbitrary GL{n, C)- 
valued matrices g(C,C), as well as for the gravitationally coupled non-linear cr-models obtained by 
dimensional reduction of Maxwell-Einstein theories in higher dimensions. For our analysis, we shall 
use the general framework of monodromy preserving deformations of ordinary differential equations 

i- 

Let us now consider the behavior of (d^'/d^)^'"^ and ((i^'/(i^)^'~^ in the complex 7-plane. Singu- 
larities in 7 arise at those points where ^'(7) is either non-holomorphic or degenerate (i.e. det ^ = 0). 
Analyticity away from the points 7 = ±1 implies that all singular points (for j = 1, N) of the 
function ^(7,^,^ regular in the sense that ^ 

^h) = G,{tmjh,C,Oil-ljf'C, as 7~7i (13) 

For 7 ^ 7j, ^'j(7,^,^) = 1 -|- 0(7 — 7j) is holomorphic and invertible. The matrices Cj and Tj are 
constant and invertible, and constant diagonal, respectively, while the (^, ^)-dependent matrices Gj 
are assumed to be invertible. The singular points depend on {(,,^) according to (^), i.e. we have 



7j = jiwj, with constants wj £ C The set {'yj,Cj,Tj} for j = 1, ...,N is generally referred 

4 



to as the set of monodromy data of ^'(7). The function ^'(7) is uniquely defined by its monodromy 
data up to normahzation |]|]. 

The logarithmic derivative ^'^^'"^ is thus holomorphic except at the points 7 = 7j where it has 
simple poles with residues 

A,{C,0 = G,T,Gf (14) 

by (^). The functions Aj{£^,^) will play a central role in the sequel. In general the number N 
of regular singularities jj may be infinite (explicit examples are the x-periodic static axisymmetric 
solutions found in |jT2|) or even continuous (this would correspond to non-constant conjugation matrices 
in the related Riemann-Hilbert problem). However, in this paper we will restrict attention to finite 
A^. Besides that, we find it convenient to impose the normalization condition ^I'^^'^^I^^oq = 0, which 
may be ensured for instance by demanding ^\^^^ = cr^. A large class of solutions with finitely many 
singularities is provided by the multisoliton solutions of Einstein's equations in ||^ (corresponding to 
matrices Tj all of whose eigenvalues are half-integer) and the finite-gap (algebro-geometric) solutions 



constructed in p3(| . 



Combining (|13| ) and (|lj) we arrive at the following differential equation in 7: 



^ = y — (15) 

57 ,tl7-7i 



Inserting (15) into (|12D, we immediately obtain 



(in the sequel summation is taken everywhere from 1 to N). Substituting (|T^ ) into (0), we get the 
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following compatibility conditions between (0) and (15) 



dAj ^ 2 ^ [Ak, A,] dA, _ 2 ^ [A,, A,\ 

e-e",^.(i-7.)(i-7,) ' dl c-e,^.(i + 7^)(i+7,) ^ ^ ^ 

It is now straightforward to check that this system is always compatible if the functions 7j obey (^). 
The first main result of this letter is thus 

Theorem 1 Let \wj G C ; j = 1,..., A^} he an arbitrary set of complex constants and Aj = 



an associated set of solutions of Then the system of linear equations ( \lq ) is always compatible, 

and the matrix function g = g{^, ^) obtained by integrating ^1^) solves 

The proof may be obtained by direct calculation. 

It is quite remarkable that the dependence of the Ernst equation and its associated linear system 
on the variables and can be completely decoupled by this theorem. In other words, the problem 
of solving Einstein's equations in this reduction has been reduced to integrating two ordinary matrix 
differential equations, which are automatically compatible unlike the original linear system (|^)! All 
information about the degrees of freedom is thereby encoded into the "initial values", i.e. the set 
of (constant) matrices A^^^ = Aj{^^^\^^^^); these are also the appropriate phase space variables, as 
we will see below. Accordingly, we will regard the functions Aj(^,^) rather than ^'(7, ^, ^) as the 
fundamental quantities from now on, and relate the system ( [T7| ) directly to the (complexified) Ernst 
equation 

Equations (O) may also be represented in "Lax form" , viz. 



dAi 



dAi 
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V\-f=^j-,Aj 



(18) 



where the matrices U and V are defined in (^ and ([l^). This form of (17) is "gauge-covariant" with 
respect to the transformation 

§ = 17(^,0*, (19) 
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Namely, the transformed function ^ satisfies the hnear system d^/d^ = U^, d^ /d^ = , where 



ii = + m^i'^ , V = n^fi'^ + nvn^^ (20) 



Clearly, the matrix functions Aj transform as Aj — > Aj = ^Aj^}^^ under (p!9[). The transformed 
matrices Aj then obey the same linear system ([l^) with the pair {U, V) replaced by {U, V). 

Theorem 1 establishes a direct correspondence between GL{2, C)-valued solutions of (|3|) and solu- 



tions of (17). However, it does not specify the conditions that must be imposed on {wj , Aj} in order 
satisfy certain restrictions which the metric g may be subject to. It is easy to see that the condition 
det 5 = 1 is guaranteed by tr^j = 0; reality of g requires the existence of an involution (complex 
conjugation) on the set {wj, Aj}. Conditions that ensure g = g"^ are more difficult to formulate and 
will be discussed elsewhere. 

Conformal factor and r-function. To each solution {Aj} of (|l^) we can associate the following 
closed one- form [||: 

9o(C,C) = E tr{AjAk)dlog{j, - 7fc) (21) 

where the exterior derivative d is to be taken with respect to the deformation parameters {(,,(,)■ The 
closure condition dqo = may be directly verified by use of ( [l7| ) and (|9|). Following the general 
prescription given in [||], we define the r-function associated with the Ernst equation by 

dlogT = qo (22) 

We will now show that this r-function has a very definite physical meaning in our context: up 
to an explicit factor, it is just the conformal factor h = e^^ [15|! To establish this result, we first 
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substitute ( |l6|) into (^); then using @ and (21) we obtain 



q = qo + 



Ij)' (1 + 7,) 



, +J2tr{AjAk)dlog{^-0 



(23) 



j<k 



Now, from ( [l7| ) it follows that J2j is {£,, ^)-independent. Furthermore, it is easy to check that 
tr^j and tr^^ are independent of ^ and ^, hence constant, for ah j. Therefore, the eigenvalues of Aj 
are (^, ^)-independent for any solution of (17) that agrees with (0). As a consequence, the expression 
J2j<k tr(^jAfc) is likewise (^, ^)-independent, and all extra terms on the r.h.s. of (p3|) may be explicitly 
integrated. 

Using (P) and (|To|), we thus arrive at 



Theorem 2 The conformal factor h ^ and the T-function (2^) are related by 



n 



dwj 



(24) 



where C & C is a constant of integration. 



Notice once more that quantities tr(^j Aj)'^ and trA| are (^, ^)-independent. If the related matrix g is 
real and symmetric, then Aj = 0, and the first factor on the r.h.s. of (24) drops out. We emphasize 
that our result is more general than previous ones (the explicit computability of h for multi-soliton 
solutions has been known for a long time [^), and valid for arbitrary non- linear <T-models coupled to 
gravity. 



Hamiltonian formulation. The system (17) is a "two-time" hamiltonian system with respect 
to the standard Lie-Poisson bracket 16 1 



{^(7) f A{fi)} = r{fi - 7) , A(7) 1 + 1® A{fi) 



(25) 



where ^(7) = ^-y\I/ ^ and the classical rational i?-matrix r(7) is equal to n/7 with 11 the permutation 
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operator in x C^. The dynamics in the and ^-directions are governed by the Hamiltonians 
H -Ok - ^ V trjA.Ak) n -Oh - ^ V tr{A,Ak) 



Compatibihty of the system (|l^) imphes {Hi^H2} = 0, i.e. the flows with respect to the two "time 
variables" ^ and ^ commute (as can also be verified by explicit computation). Note that our formulation 
is far simpler both technically and conceptually than previous hamiltonian treatments of such systems, 
which were based on the use of "one-time" Hamiltonians, and where the Lie-Poisson brackets (|25| ) 
would in addition depend on the space coordinates. The notorious problems caused by derivatives of 
(5-functions (so-called "non-ultralocal" terms) in the relevant Poisson brackets are altogether avoided 
here. Furthermore, the cumbersome structure of the canonical current algebra in the conventional 
approach is replaced by a more transparent algebraic structure in our framework. 

To be sure, we should regard (^6|) as constraints a la Dirac rather than conventional Hamiltonians, 
because (|3|) is derived from a generally covariant theory. To do this properly would, however, require 
that we undo the choice of Weyl coordinates, on which ([l|) and (|3|) are based, and to treat ^ and ^ 
as canonical variables subject to {2A;^,^} = {2k^,^} = 1 and {2k^,^} = {2/c|-, ^} = 0. The quantities 
$1 := 2k^ — Hi and <I>2 := — H2 are thereby converted into (mutually commuting) constraint 
operators on an enlarged phase space (they are, in fact, just the Virasoro constraints). The "time 
evolutions" of the spectral parameter 7 are also generated canonically in the sense that 7^ = {<I>i,7} 
and 7^ = {^*2)7} (so "time" must eventually be quantized in this scheme!). 

Final remarks. 

1. An obvious advantage of using the variables Aj in comparison with the ones employed tra- 
ditionally in this context is that they generate a closed Lie algebra with respect to the standard 
Lie-Poisson bracket (25). Secondly, the common features of our system ( |T7|) with the classical limit 



of the Knizhnik-Zamolodchikov equations 17| suggest that one should quantize (|T^ in analogy 
with the KZ equations (although quantum gravity will certainly introduce new features). We have 



only discussed the case of finite in this paper, but there are in principle no obstacles to considering 
N = oo from the outset, since the space of finite N solutions can be naturally embedded in this 
larger space. For illustrative purposes, it is quite useful to think of the set of finite solutions as 
the 'W-particle sector" of the theory because, depending on the reality conditions, the solutions g 
corresponding to ([l5|) possess exactly N or singularities wi,...,wn in the upper ^-half-plane. In 
fact, a proper treatment of the Ernst equation as a (generally covariant) quantum field theory will 
presumably require taking into account as a "particle number operator" . 

2. The extension of our results to the case of a Lorentzian world-sheet and to arbitrary G/H coset 
space cj-models coupled to gravity in two space-time dimensions is straightforward. In the notation 
of where this case is reviewed in some detail, the linear system matrix ^' corresponds to Vr/(V)^^, 
where 77 denotes the Cartan involution (e.g. r/(V) = (V"^)~^ for G = SL[n)); furthermore, the spectral 
parameter t used there corresponds to i times the parameter 7 employed in the present paper. For 
arbitrary coset spaces, the matrices Aj belong to the Lie-algebra of G; like g, they may be subject to 
further restrictions. 

3. Analogs of the static axisymmetric (multi-Schwarzschild) solutions for arbitrary cj-models can 
be easily constructed in our formalism by choosing the matrices Aj in the Cartan subalgebra of the 
relevant Lie algebra. From (|l^) it is then immediately evident that Aj = const. 



4. Obviously our formulation will yield a new realization of the Geroch group [18| and its gener- 
alizations; we here just note that this group mixes sectors belonging to different "particle numbers". 
It is known that the corresponding Kac Moody algebras act on the conformal factor via their central 
extension |1£, ^, |9|; combining this result with our Theorem 2 should shed some light on the group 
theoretical meaning of the r-function. 

A detailed account of the results described in this letter is in preparation. 
Acknowledgements. We thank Alexander Kitaev for important discussions. 



10 



References 

[I] M. Jimbo, T. Miwa and K. Ueno, Physica ID (1981) 80 

[2] A. C. Newell, Solitons in Mathematics and Physics, Philadelphia, SIAM (1985) 

[3] L. Takhtajan and L. Faddeev, Hamiltonian Methods in the Theory of Solitons, Springer Verlag 
(1986) 

[4] V. Korepin H. Bogolubov and A. Izergin, Quantum Inverse Scattering Method and Correlation 
Functions, Cambridge Univ. Press (1993) 

[5] D. Kramer, H. Stephani, E. Herlt and M. MacCallum, Exact Solutions of Einstein's Field equations, 
Cambridge University Press (1980) 

[6] D. Maison, Phys.Rev.Lett. 41 (1978) 521 

[7] V. A. Behnskii and V. E .Zakharov, Sov.Phys. JETP 48 (1978) 985 
[8] P. Breitenlohner and D. Maison, Ann. Inst. Poincare 46 (1987) 215 

[9] H. Nicolai, in Recent Aspects of Quantum Fields, eds. H. Mitter and H. Gausterer, Springer Verlag 
1991. 

[10] Analogous expressions were obtained in: 

G. Neugebauer and D. Kramer, J.Math.Phys. 32 (1991) 2727 

[II] By contrast, integrable systems with constant spectral parameter such as the non-linear 
Schrodinger and Sine-Gordon equations will always exhibit essential singularities in the 7-plane. 
Similarly, the crucial relations ( |T2|) have no analog for systems with constant spectral parameter. 

[12] R. C. Myers, Phys.Rev. D35 (1987) 455; 

D. Korotkin and H. Nicolai, Nucl.Phys.B429 (1994) 229 



11 



[13] D. Korotkin, Theor.Math.Phys, 77 (1989) 1018; 

D. Korotkin and V. Matveev, St.Petersburg Math. J. 1 (1990) 379 

[14] Treating as independent parameters, we would get the so-caUed Schlesinger equations instead 
of ([l|), see: 

L. Schlesinger, J.Reine u. Angew. Math. 141 (1912) 76. 

[15] The close link between r and h was conjectured some time ago by A.Kitaev (private communi- 
cation). 

[16] E. Sklyanin, J.Sov.Math. 47 (1989) 2473; 

J. Harnad, Quantum isomonodromic deformations and the Knizhnik-Zamolodchikov equations, 



preprint hep-th/940607S, 1994 



[17] N. Reshetikhin, Lett.Math.Phys. 26 (1992) 167 
[18] R. Geroch, J.Math.Phys. 12(1971)918 

[19] B. Julia, in Unified Theories and Beyond, Proc. 5th Johns Hopkins Workshop on Current Prob- 
lems in Particle Theory, Johns Hopkins University, Baltimore (1980) 



12 



